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O ■ Abstract. 

S ! Comprehensive set of shell model calculations for ^^"^''Se isotopes have been 

performed with recently derived interactions, namely JUN45 and jj44b for /5/2P59/2 
space. To study the importance of the proton excitations across Z — 2& shell in this 
^ I region mentioned by Cheal et al. [Phys. Rev. Lett. 104, 252502 (2010)], calculation 

^»0 ' for fpgg/2 valence space using an fpg effective interaction with ^^Ca as core and 

^jD ^ imposing a truncation has also been performed. Comparison of the calculations with 

experimental data show that the predicted results of jj44b interaction are in good 
agreement with experimental data. 



PACS numbers: 21.60. Cs, 27.50. +e 



1. Introduction 

The nuclei near Z = 28 region are subject of intensive research of both experimental 
and theoretical investigations IH |2l |3]. In particular, the region of neutron- rich nuclei 
around A^ = 40 to A^ = 50 shows an evolution of shell structure. The nucleon-nucleon 
interaction, the spin-orbit, tensor part, and three-body effect play an important roles 
in the shell evolutions. Due to tensor interactions, the nuclear mean field undergoes 
variations with neutron excess. This leads to monopole migration, which is observed 
in both side of stability line. As we approach towards the neutron drip line, the 
neutron density becomes very diffused, which can also leads to shell quenching. Further 
motivations for nuclear structure studies are for understanding the other open questions 
like, how does the filling of neutron orbitals beyond N = 28 around ^^Ca influence the 
shell structure, what is the role of tensor part of proton-neutron interaction, how does 
the spin-orbit splitting between f^/2 and /7/2 evolve when approaching A^ = 50, and 
what is the specific role of particle-hole excitations through Z = 28 and A^ = 40 shell 
gaps for the onset of deformation below and above ^^Ni. 
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Figure 1. Systematic of the experimentally observed E{2f) for Z — 28 to Z — 34 
near the N — 28, 40 and 50 shell closure. 



To study the proton excitations across Z = 28 gap many experimental 
investigations have been performed for Cu [4J , Ga [5J and As [6j isotopes. The inversion 
of 77/5/2 and 77^3/2 orbitals for Cu isotopes have been established while measuring 
magnetic moments for the ground-state. In case of ^^~®^Ga isotopes experimental 
measurement of spin and moments reveal that there is structure change between A^ = 40 
and A^ = 50. The 77/5/2 orbital dominant in the g.s. of ^^Ga and for ^^Ga, 5/2" level 
become ground state. The evolution of structure can be seen from emptying of 77^3/2 



orbital to t^J^j^ is started as we move from ^^Ga to '''^Ga. Finally F=5/2~ become 
ground state for *^Ga. In case of ™Ga, the f^/2P99/2 space fail to explain correctly the 
magnetic and electric quadrupole moment. Thus, it is crucial to include 7r/7/2 orbital 
in the model space to study the effect of proton excitations across the Z = 28 effect. 
In case of As {Z = 33), recent experimental investigation by Astier et al [6] suggests 
that f5/2Pg9/2 space is not enough to explain quadrupole excitation built on the 5/2^ 
and 9/2t state of ^^As. 

In this work, shell model calculation for ^^~*^Se isotopes in two different model 
spaces f5/2Pg9/2 and fpg9/2 have been performed. The aim of this study to study the 
importance of proton excitations across Z = 28 gap. One of us recently reported 
importance of 7rfi/2 orbital for predicting moments of Ga isotopes [7]. In Fig. [H the 
systematics of 2"'" for Ni to Se are shown. This figure shows the persistence of A^ = 50 
shell closure while A = 40 disappears as we move from Ni to Se. 

In section [2] the details about shell model calculations is described and then 
spectroscopic results for even-even 78,80,82,84gg ^^^^ presented in section [3l In section 
4, transition probability, quadrupole moments and occupation numbers are presented. 
Finally section |5] gives the concluding remarks. 

2. Details of model spaces and interactions 

The present shell model calculations have been carried out in the f 5/2 P 99/2 and f pg9/2 
spaces. In the /5/2P 5*9/2 valence space the calculations have been performed with the 
interactions JUN45 [8] and jj44b [9]. The single-particle energies for the 1^3/2, O/5/2, 
lpi/2 and 0^^9/2 single-particle orbits employed in conjunction with the JUN45 interaction 
are -9.8280, -8.7087, -7.8388, and -6.2617 MeV respectively. In the case of the jj44b 
interaction they are -9.6566, -9.2859, -8.2695, and -5.8944 MeV, respectively. The core 
is ^^Ni, i.e. N = Z = 28, and the calculations are performed in this valence space 
without truncation. 

In the f pg9/2 valence space, we use a ^®Ca core, i.e. only the protons are active 
in the fi/2 orbital, and the interaction fpg reported by Sorlin et al [10]. The single- 
particle energies are 0.0, 2.0, 4.0, 6.5 and 9.0 MeV for the O/7/2, 1^3/2, lpi/25 O/5/2, and 
0(^9/2 orbits, respectively. Since the dimensionality of this valence space is prohibitively 
large, we have introduced a truncation by allowing t^r particle-hole excitations from the 
7r/7/2 orbital to the upper fp orbitals (^3/2, /5/2 5 P1/2) for protons and ty particle-hole 
excitations from the upper fp orbitals to the i'g9/2 orbital for neutrons. The maximum 
allowed value for t^ and ty is four. The maximal dimension 165000813 ( ~ 10^ ) is 
reached for positive parity in the case of ^^Se when using f5/2pg9/2 space with ^^Ni core 
since neutron number is furthest from the closed shell for this nucleus among the Se 
isotopes considered in this work. In case of ^^Se computing time ~ 21 days for both 
parity. The calculations were performed with shell-model code ANTOINE pT] . 



3. Spectra analysis 

The results for the three interactions for different model spaces used in the calculations 
are presented with respect to the experiment in Figs. 2-5. Earlier shell model result in 28- 
50 model space for pairing plus quadrupole-quadrupole interaction has been reported 
in literature by Yoshinaga et al. [12]. Further this work will add more information 
to the earher work p[2], by including /7/2 orbital in the model space to study the 
proton excitation across Z = 28 shell. The shell model structure of Se isotopes was 
discussed in [8] using JUN45 interaction. The results for the three interactions used in 
the calculations are presented with respect to the experiment. 

3.1. ^^Se 

Comparison of the calculated values of the energy levels of ^^Se with the experimental 
data is shown in Fig. |2l All the three interactions correctly reproduces the experimental 
sequence of 0^, 2^, 4^, 6^, 8^, 10^ and 12^ levels at 614, 1503, 2546, 3585, 4625 
and 8784 keV, respectively. In case of JUN45, the calculated 25*", Ai and 61" levels are 
134, 240 and 251 keV higher than the experimental ones, while remained 8^, lO]*" and 
12 j'" are predicted 117, 323 and 397 keV lower than experimental levels. The 25*" level 
calculated by jj44b is only 60 keV higher as compared to its experimental counterpart. 
Difference between experimental and calculated energies is increased by increasing the 
excited energies reaching 433 keV at 12^. All the levels calculated by fpg interaction 
is much lower than in the experiment. By the increasing spin the difference between 
the calculated and experimental levels is increased and it become rather 1751 keV when 
reaching 12^ level. 

In all the calculations first two 0^ and 2+ levels of the second positive-parity band 
are interchanged with respect to experimental ones. The next O"*" and 3"'" levels are 
located higher than in the experiment in all calculations. The experimental sequence 
levels 2+, 4+, 4+ and 5+ at 1996, 2191, 2682 and 2735 keV are better reproduced by the 
JUN45 calculation. The two 6^ levels appear in all calculations after the 5"*" level which 
have not been observed in the experiment. The last three levels of this band calculated 
by the JUN45 interaction excellently agree with the experimental data. Experimental 
levels of this band measured up to 10"*" and many odd spin levels 1"*", 7"*", 9"*", which 
appear in the calculations, have not been observed in the experiment. 

The experimental negative-parity levels start with the 3~ level. The first negative- 
parity level is 4~ in both JUN45 and jj44b calculations, while it is 2~ in the fpg 
calculation. In the JUN45 calculation the first few levels are much compressed than 
in the experiment. The spins and parities of (7^), (8^), (9~), and (9~) levels which 
is not confirmed in the experiment, are predicted by all calculations as negative parity 
levels. The sequence of experimentally observed levels starting from 1~ is exactly the 
same as in the jj44b calculation. 



3.2. «0Se 

In Fig. [3] we have shown the comparison of the values of the energy levels calculated by 
JUN45, jj44b and fpg interactions with experimental data. The values calculated by 
JUN45 are in very good agreement with the all levels of the first experimental positive- 
parity band which are available up to 10"*". The jj44b calculation predicts higher values 
after 4"*", while the values predicted by the fpg calculation are lower up to 4"*" and then 
the G"*" level is higher. The 8^ and lO"*" levels are again lower than in the experiment. 

The fist two levels of the second positive-parity band are interchanged with respect 
to those of ^^Se in jj44b calculation. The sequence of these levels are not changed 
with respect to those of ^^Se in JUN45 calculation. In the fpg calculation 2"*" level is 
lowest in the second positive-parity band. The 4"*" levels is predicted much lower than 
in the experiment, while the second O"*" and 2"*" levels are higher in this band than in 
the experiment. The second 4"*" level is only 17 keV higher than in the experiment. 
Then three G"*" levels in succession come in the calculation. The S"*" in the calculation 
is much higher than in the experiment. For odd spin-positive-parities, which we have 
shown in the third column, all the calculation predict two levels for each spin (l"*" and 
l"*", S"*" and S"*",..). In the experiment only one l"*" and S"*" are measured, first of which is 
higher than all l"*" levels in all calculation, while second of them lower than S^" level in 
all calculations. 

The lowest negative-parity level according to JUN45 and jj44b calculations are 1~, 
which is much higher in the experiment, however the second and the third 1~ levels are 
closer to the experimental ones in both calculations. Four 3~ levels are measured in 
the experiment. For both JUN45 and jj44b calculations the reported three 3~ levels are 
located approximately the same in both calculations. We have also reported three 5~ 
levels in all calculations. All of them are lower than experimentally measured two 5~ 
levels. 

3.3. 82 Se 

As is seen from Figure H] in all calculated levels in the first column are in similar pattern 
to experimental one. Better values are predicted by JUN45 calculation. 

In all calculations the first level in the second column is O"*" level like in the 
experiment. Closest value of this level is predicted by JUN45. The spacing between 
the first and the second 2"*" levels in the second column are less than that of in the 
experiment in all calculations. In JUN45 calculation first of two 4+ levels is located 
lower while the second one is higher than in the experiment. In jj44b calculation both 
of them are much higher than in the experiment. In fpg calculation the first 4"*" is only 8 
keV higher than in the experiment. The second 0"^ is predicted well by jj44b calculation. 

Among the positive-parity odd spin levels only 9^ and (11) level energies measured 
in the experiment. The (11) level's parity is predicted to be positive by the calculations. 
As in the case of previous isotopes there are many other 1"*", 3+, 5"*" levels in the 
calculation which are not measured in the experiment. 
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As is seen from the Figure H] agreement of the all calculated values of negative- 
parity levels with the experimental data is improved very much as compared to the 
78,80gg isotopes. 

3.4. 84Se 

The levels of the first column of Figure [5] are predicted better by jj44b calculation. Also 
G"*" level is little bit higher than the in JUN45, spacings between the levels very much 
like to the experimental ones. 

In all calculations the first positive-parity level in the second column is O"*" as in the 
experiment. The sequences of Of, 02" and 2f, 23" levels is the same with the experimental 
one in jj44b calculation and spacings between them are larger than in the experiment. 
Sequence of the pair of these levels are different from the experimental one in JUN45 
and fpg calculations and spacing between them are much larger than in the experiment. 
The 4f , and 4^ levels are located much lower than in the experiment in all calculations. 

The l"*" and 3"^ levels appearing in the calculations have not been measured in the 
experiment. The measured S"*" level is predicted better by JUN45. The 7~^ level is 
located very high in all the calculations as compared to experimental one. 

The structure of the negative-parity levels for ^^Se is changed so drastically as 
compared to previous isotopes and both JUN45 and jj44b calculations fail to explain it. 



4. Transition probability, quadrupole moments and occupation numbers 
analysis 

4.1. E2 transition probability and quadrupole moments 

The calculated B{E2) transitions are shown in table 1. In case of ^^'^°Se, the results 
predicted by fpg interaction show better agrement with experimental data. While for 
^^'®^Se, the jj44b interaction results are more reasonable. Thus we may conclude that 
proton excitation across Z = 28 shell for lighter Se isotopes are important. We have also 
calculated static quadrupole moments as shown in table 2. For the first 2+ state, the 
JUN45 and jj44b interaction correctly predict sign of quadrupole moments for 78,80,82gg_ 
The fpg interaction predicts positive sign for ^^Se. The overall agreement for quadrupole 
moments predicted by jj44b is better. 

4-2. Occupation numbers 

In Fig. 6, we show the proton/neutron occupation numbers of fpg-shell orbits for Of 
and 2f levels. The proton occupancies increasing smoothly in the /5/2 orbital, while P3/2 
decreasing. But beyond A^ = 46, in case of JUN45 interaction the change in occupancy of 
these two orbitals are very significant. With both interaction the occupancy of proton 
Pi/2 and gg/2 orbitals are similar. As the neutron number increases the occupation 
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Figure 2. Comparison of experimental and calculated excitation spectra of ^^Se with 
three different interactions. 
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Figure 3. Comparison of experimental and calculated excitation spectra of ^''Se with 
three different interactions. 
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Figure 4. Comparison of experimental [13] and calculated excitation spectra of ^^Se 
with three different interactions. 
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Figure 5. Comparison of experimental [T4l[T5l[T6] and calculated excitation spectra 
of ^^Se with three different interactions. 



Table 1. B{E2) reduced transition strength in W.u. Effective charges Cp 
e„ — 0.5 were used. Experimental values were taken from the NNDC database. 
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= 1.5 





78 Se 


80Se 


82Se 


84Se 


BE(2+ ^ 0^) 










Experiment 


32.8(4.5) 


24.2 (0.4) 


16.7 (0.3) 


N/A 


JUN45 


18.83 


16.53 


13.88 


6.64 


jj44b 


20.74 


18.88 


15.33 


8.37 


fpg 


30.87 


24.08 


16.33 


8.22 


BE(4+ -> 2^) 










Experiment 


38.2t|:? 


34.7 (1.1) 


18.7 (3.0) 


N/A 


JUN45 


25.47 


21.40 


19.60 


1.96 


jj44b 


27.23 


25.65 


20.77 


0.005 


fpg 


46.27 


32.85 


25.86 


1.42 


BE(6+ ^ 4+) 










Experiment 


48(14) 


N/A 


N/A 


N/A 


JUN45 


23.69 


16.81 


15.42 


1.02 


jj44b 


24.22 


21.92 


18.01 


3.25 


fpg 


48.92 


32.01 


24.42 


0.12 


BE{8t -^ Qt) 










Experiment 


57(19) 


N/A 


0.56 (0.03) 


N/A 


JUN45 


13.01 


7.45 


0.38 


0.0008 


jj44b 


21.86 


1.27 


0.31 


0.002 


fpg 


49.01 


7.13 


0.57 


1.13 



number of iygg/2 orbital increases drastically. The neutron number occupation show a 
similar distribution for 2f as in the ground state. 
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Table 2. Electric quadrupole moments, Qs (in eb), with the three different 

interactions (the effective charges ep=1.5, e„=0.5 are used in the calculation). 





^^Se 


soSe 


«2Se 


s^Se 


Q(2|) 










Experiment 


-0.20 (7) 


-0.31 (7) 


-0.22 (7) 


N/A 


JUN45 


-0.13 


-0.31 


-0.33 


+0.01 


jj44b 


-0.32 


-0.36 


-0.37 


-0.27 


fpg 


+0.47 


-0.35 


-0.36 


+0.04 


Q(2+) 










Experiment 


+0.17 (9) 


N/A 


N/A 


N/A 


JUN45 


+0.13 


+0.28 


+0.24 


-0.09 


jj44b 


+0.30 


+0.35 


+0.27 


+0.13 


fpg 


-0.33 


+0.36 


+0.28 


+0.007 


Q(4|) 










Experiment 


-0.68 (15) 


N/A 


N/A 


N/A 


JUN45 


-0.09 


-0.35 


-0.39 


+0.10 


jj44b 


-0.36 


-0.40 


-0.42 


+0.18 


fpg 


+0.63 


-0.29 


-0.43 


+0.15 



5. Summary 

In summary, comprehensive study for the structure of neutron-rich even-even Se iso- 
topes have been carried out using large-scale shell-model calculations for two spaces: 
full /5/2Pfi'9/2 space and fpg9/2 space with ^®Ca core. 



The following broad conclusions can be drawn: 

The overall calculated results for the energy levels, B{E2)s and quadrupole 
moments are in good agreement with the experimental data. 

The E2 transitions, quadrupole moments analysis show the importance of proton 
excitations across Z = 28 shell for fpgg/2 space. 

For ^^~^°Se the B[E2) values predicted by fpg transitions are in better agreement 
with experimental data. 

The result (wave functions) of ^^Se, may be used for calculating nuclear transition 
matrix elements and finally half-live for this nucleus which is a good candidate for 
neutrinoless double beta decay. 

It is also important while tuning the effective interaction we can also take 
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Figure 6. (Color online) Proton/Neutron occupation numbers of the JUN45 and 
jj44b (p3/2, /5/2, Pi/2 and ^9/3 -shell orbits) and jpg {hji-, P3/2, /5/2, P1/2 -shell 
orbits) interactions- for two low-lying states in even-even Se isotopes. (Upper panel) 
Oi states; (lower panel) 2f states. 



experimentally known quadrupole moment as a parameter to increase predictive 
power of effective interaction. 

Further theoretical development is needed by enlarging model space by including 
z/(i5/2 orbital to study simultaneously proton and neutron excitations across Z = 28 
and Z = 50 shell, respectively. 
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